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The current important issue in numerical relativity is to determine which formulation of the 
Einstein equations provides us with stable and accurate simulations. Based on our previous work on 
"asymptotically constrained" systems, we here present constraint propagation equations and their 
eigenvalues for the Arnowitt-Deser-Misner (ADM) evolution equations with additional constraint 
terms (adjusted terms) on the right hand side. We conjecture that the system is robust against vio- 
lation of constraints if the amplification factors (eigenvalues of Fourier-component of the constraint 
, propagation equations) are negative or pure-imaginary. We show such a system can be obtained 

■ by choosing multipliers of adjusted terms. Our discussion covers Detweiler's proposal (1987) and 

^S) ' Frittelli's analysis (1997), and we also mention the so-called conformal-traceless ADM systems. 

! PACS numbers: 04.20.-q, 04.20.Fy, 04.25.-g, 04.25.Dm 



I. INTRODUCTION 

\ The effort to solve the Einstein equations numerically - so-called Numerical Relativity - is now providing an inter- 
■ esting bridge between mathematical relativists and numerical relativists. Most of the simulations have been performed 
' using the Arnowitt-Deser-Misner (ADM) formulation or a modified version. However, the ADM formulation has 
, not been proven to be a well-posed system, since its evolution equations do not present a hyperbolic form in its 

original/standard formulation. 
^ Most simulations are performed using the "free evolution" procedures: (1) solve the Hamiltonian and momentum 
constraints to prepare the initial data, (2) integrate the evolution equations by fixing gauge conditions, and (3) monitor 
the accuracy /stability by evaluating the constraints. Many trials have been made in the last few decades, but we 
O . have not yet obtained a perfect recipe for long-term stable evolution of the Einstein equations. Here we consider the 
^ '' problem through the form of the equations. 

5-H One direction in the community is to rewrite the Einstein evolution equations into a hyperbolic form and to apply it 

to numerical simulations |Q . This is motivated by the fact that we can prove well-posedness for the evolution of several 
J> , systems if they have a certain kind of hyperbolic feature. The authors recently derived ||,^ three levels of hyperbolic 
system of the Einstein equations using Ashtekar's connection variables and compared them numerically 

. , We found that (a) the three levels of hyperbolicity can be obtained by adding constraint terms and/or imposing 
' gauge conditions, (b) there is no drastic difference in the accuracy of numerical evolutions in these three, and (c) the 
symmetric hyperbolic system is not always the best for reducing numerical errors. Similar results regarding to (a) 
and (b) are reported by Hern based on the Frittelli-Reula formulation § . 

What are, then, the criteria for predicting the stable evolutions of a system? Inspired by the "A-system" proposal 
we have considered a so-called "asymptotically constrained" system, that is, a system robust against the violation 
of the constraints [|lO|. The fundamental idea of the "A-system" is to introduce artificial flow onto the constraint 
surface. However, we also found that such a feature can be obtained simply by adding constraint terms to the 



evolution equations which we named "adjusted systems" |11|. We explained the reason why this works by analyzing 



the evolution equations of the constraints (the propagation of the constraints) . We proposed that the stablity of the 
system can be predicted by analyzing the eigenvalues (amplification factors) of the constraint propagation equations 
(We describe this in detail in §|l|) . We confirmed that our proposal works both in Maxwell and Ashtekar systems . 



'^We derived weakly, strongly (= diagonalizable) and symmetric hyperbolic systems. The mathematical inclusion relation is 

weakly hyperbolic 9 strongly hyperbolic 9 symmetric hyperbolic. 

See details in Ml. 
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The purpose of this article is to apply our proposal to the ADM systeni(s). Especially, we consider the "adjusting 
process" (adding constraints in RHS of evolution equations), and the resultant changes to the eigenvalues of the 
constraint propagation systems. This adjusting process can be seen in many constructions of hyperbolic systems in 
the references. In fact, the standard ADM for numerical relativists is the version which was introduced by York 
where the original ADM system has already been adjusted using the Hamiltonian constraint (see more detail in 
§111). The advantage of the standard ADM system is reported by Frittelli [|l3| from the point of the hyperbolicity 
and the characteristic propagation speed of the constraints. Our discussion extends her analysis to the amplification 
factors. 

One early effort of the adjusting mechanism was presented by Detweiler Our study also includes his system. 



jtuay 

and shows that this system actually works as desired for a certain choice of para meter (§|rv|). We also study the same 



procedure for the "conformal-traceless" ADM (CT-ADM) formulations |15 1^ which is recently the most popular 
system in numerical simulations (§[v|). 

The analysis in the text is for perturbational violation on a flat background. Further applications are available, but 
we will discuss them in future reports. In the Appendix, we also give numerical demonstrations of the adjusted- ADM 
systems discussed in the text. 



II. CONSTRAINT PROPAGATION AND "ADJUSTED SYSTEM" 



We begin by reviewing the background of "adjusted systems" and our conjecture. 

The notion of the evolution equations of the constraints is often discussed from the point of whether they form a 
first class system or not. Fortunately, the constraints in the (original/standard) ADM formulation arc known to form 
a first class system. Due to this fact, numerical relativists only need to monitor the violation of the Hamiltonian and 
momentum constraints during the free evolution of the initial data. 

Our essential idea here is to feed this procedure back into the evolution equations. That is, we adjust the system's 
evolution equations by characterizing the constraint propagation in advance. Let us describe the procedure in a 
general form. Suppose we have a set of dynamical variables, M°(a;*,t), and its evolution equations, 

ft^^" = /(^i^9,^i^•••), (2.1) 

which should satisfy a set of constraints, C''(u°, diu'^, • • •) w 0. The evolution equation for C° can be written as 

dtC = g{CP,d,CP,---). (2.2) 



We can perform two main types of analysis analysis on (2.2) 



1. If (2.2) is in a first order form (that is, only includes first-order spatial derivatives), then the level of hyperbolicity 
and the characteristic speeds (eigenvalues A' of the principal matrix) will definitely deterine the stability of the 
system. We expect mathematically rigorous well-posed features for strongly or symmetric hyperbolic systems, 
and the characteristic speeds suggest to us satisfactory criteria for stable evolutions if they are real, and under 
the propagation speed of the original variables, m", and/or within the causal region of the numerical integration 
scheme applied. 



2. On the other hand, the Fourier transformed- (2. 2) 



dtC''=g{CP), (2.3) 

where CP{x,t) = JCP{k,t) exp{ik ■ x)d^k, also characterizes the evolution of the constraints independently of 
its hyp erbolicity. As we have proposed and confirmed in jl^ , the set of eigenvalues A' of the coefficient matrix 



in (2.3) provides a kind of amplification factor of the constraint propagation, and predicts the increase/decrease 
of the violation of the constraints if it exists. More precisely, we showed in that 



if the eigenvalues of (2.3) (a) have a negative real-part, or (b) are non-zero (pure-imaginary) eigenvalues, 



then we see more stable evolutions than a system which does not. 

This is because the constraints are damped if the eigenvalues are negative, and are propagating away if the 
eigenvalues are pure imaginary. We found heuristically that the system becomes more stable (accurate) when as 
much As satisfies the above criteria and/or as large magnitude of As away from zeros. (Examples in pj] ] are of 
the plane wave propagation in the Maxwell system and the Ashtekar system.) We remark that this eigenvalue 
analysis requires that we fix a particular background metric for the situation we consider, since the amplification 
factor depends on the dynamical variables 
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The above features of the cons trai nt propagation, (2.2), will change when we modify the original evolution equations. 
Suppose we adjust the RHS of {2A) by adding the constraints. 



dtu^ = f{u-, d,u\ ■■■)+ FiC, d,cp, ■ ■ •), 



then (Ej2) will also be modified as 



dtCP = giC, d^C, • • •) + G(C^ d.C, ■ ■ •). 



(2.4) 



(2.5) 



By taking the characteristic speed of ( p.5[ ) and the amplification factor of the Four ier transformed- ( p. 5[ ), the predicted 
stability of the system (2.4) becomes different to that of the original system, (2.2). 



Our proposed "adjusted system" is obtained by finding a certain functional form of F(C, diC , ■ ■ ■) in (2.4) so as 
to get a more stable prediction in the analysis of the eigenvalues A' and A*. In the following discussion, we show two 
eigenvalues A' and A ' for each ADM system. We remark again that the term 'characteristic speed' here is not for the 
dynamical equation (|2.1|), but for the constraint propagation equations (2.2). 



III. STANDARD ADM SYSTEM 



A. Standard ADM system and its constraint propagation 



We start by analyzing the standard ADM system. By "standard ADM" we mean here the most widely adopted 
system, due to York [n2| , with evolution equations 



(3) 



and constraint equations 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



where {jij,Kij) are the induced three-metric and the extrinsic curvature, {a, Pi) are the lapse function and the shift 
covector, is the covariant derivative adapted to 7^^, and i?-^- is the three- Ricci tensor. 

The constraint propagation e quat ions, which are the time evolution equations of the Hamiltonian constraint (|3.3| ) 
and the momentum constraints (3.4), can be written as 



dtH = P'{d,n) ~ 2af'id,Mj) + 2aKn + a{dnmk){2j"''-f''' - 7"S'')-Mj - if^^Ja)M^, (3.5) 
dtM, = -il/2)a{d,n) + P^djMi) + aKM, - {d,a)n - (3"-/^' {daikjM, + {dAh^'Mj. (3.6) 

The simplest derivation of ( |3.5|) and ( |3.6[ ) is by using the Bianchi identity, which can be seen in Frittelli (Note 
that C in is half our 7i, and we have corrected typos in eq.(ll) in [p3[). 
The characteristic part of ( [3.5| ) and (3^) can be extracted as 



(n 



-{l/2)aSl P'Sf 



2af'\ f n 



n 



\M,) - '\M, 

which indicates that the characteristic speeds (eigenvalues of the characteristic matrix, P') are 

A' = (/?',/?', ± a\/7") (no sum over 0- 

Since rank(P' — /3') = 2, the matrix P' is diagonalizable, but not the symmetric. 
Simply by inserting (1/2) in front of Ti. above, we obtain 



('HI2\ 



— a7 
P'5l 



(3.7) 



(3.8) 



(3.9) 



the characteristic matrix becomes symmetric (with the same eigenvalues). This is a feature of the standard ADM 
system that was pointed out by Frittelli. (Actually 7i/2 is the form originally given by the Lagrangian formulation.) 
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B. Amplification factors on the Minkowskii background 



As a first example, we consi der the pe rtur bation of the Minkowskii spacetime: a — 1, (3^ = 0, jij — Sij. By taking 
the Hnear order contribution, ( |3.5| ) and (3.6) are reduced to 



dt 



-2ikj 
-{l/2)ih 



(3.10) 



in Fourier components. The eigenvalues of the coefficient matrix of (3.1C), which we call amplification factors, become 

A' = (0,0,±iVfc2), (3.11) 

where fc^ = k'^ + ky + k^. These factors will be compared with others later, but we note that the real parts of all the 
As are zero. 



IV. ADJUSTED ADM SYSTEMS 
A. Adjustments 



Generally, we can write the adjustment terms to (3.1) and ( |3.2[ ) using (3.3) and (3.4) by the following combinations 
(using up to the first derivative of constraints), 

adjustment term of : +Pij'H + Q'^ ijMk + p'' ijiDkU) + q''\j{DkMi), (4.1) 

adjustment term of dtK.j : ^R^jU + S^^^Mk + r^r, {DkU) + s^\j {DkMi), (4.2) 

where P, Q, R, S and p, q, r, s are multipliers (please do not confuse Rij with three Ricci curvature that we write as 
Rij )■ Since this expression is too general, we mention some restricted cases below. 

We remark that our starting system, ( pTl| ) and ( p.2[ ), is the standard ^DM system for numerical relativists introduced 
by York |]l^ . This expression can be obtained from the originally formulated canonical expression by ADM jl] , but in 
that process the Hamiltonian constraint equation is used to eliiminate the three dimensional Ricci scalar. Therefore 
the standard ADM is already adjusted from the original ylDM system. We start our comparison with this point. 



B. Original ADM vs Standard ADM 



Frittelli's adjustment analysis ||13| can be written in terms of ( |4.l| ) and ([4.2|), as 

% = (1/4)(a'- l)a7^J, 



(4.3) 



where /Lt is a constant and set other multiplier zero. Here /i = 1 corresponds to the standard ADM (no adjustment, 
since Rij = 0), and fj, = to t he original ADM (without any adjustment to the canonical formulation by ADM). 

Keeping the multiplier (4.3) in mind, we here discuss the case of non-zero Rij, S^ij (and all other multipliers zero) 
case. The constraint propagation equations become 



dtU = P'{d,n) - 2af={d,Mj) + 2aKH + a(9nmfe)(27"'7'^' - i"'''i')Mj - A^'^{dja)M, 
+2KRn - 2K'^Rij-H + 2K-i'^S^\jMk - 2K'^S''^jMk, 

-~{l/2)a(d,H) + f3^{djM.,) + aKM^ - {d,a)n - (3''-/^' {d,-fik)Mj + {OAh^'Mj 
+^''^{^,Rk^)■H - ^'\d,R,k)n + Rh{d,n) - Rjkl'\d^n) 



(4.4) 



dtM 



that is, (4.4) and (|4.5|) form a first-order system. The principal part can be written as 



(4.5) 
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n 



(4.6) 



The general discussion of the hyperboUcity and characteristic speed of the system (4.6) is hard, so hereafter we restrict 
ourselves to the case 



where we recover (4.3) by choosing ki = (/i — l)/4 and K2 = 0. The eigenvalues of (4.6) then become 
A' = (1 - K2)/3' ± y^a27«(l + 4ki) + (^2/3')^) (no sum over I) 



(4.7) 



(4.8) 



and the hyperbolicity of (4.6) can be classified as (i) symmetric hyperbolic when ki = 3/2 and K2 — 0, (ii) strongly 
hyperbolic when a'^'j^''{l + Aki) + k|(/3')2 > where ki 7^ —1/4, and (iii) weakly hyperbolic when a^'j'-''{l + Aki) + 

4(0'? > 0. 

For the case of (4.7) on a Minkowskii background metric, the linear order terms of the constraint propagation 
equations become 



d. 







^ 



-(l/2)(l + 4Kl)^fc, 
whose Fourier transform gives the eigenvalues 

A' = (0,0,±V-A:2(1+4ki)) 



(4.9) 



(4.10) 



That is (two Os, two pure imaginary) for the standard ADM, and (four Os) for the original ADM system. Therefore, 
according to our conjecture, the standard ADM system is expected to have better stability than the original ADM 
system. 



C. Detweiler's system 

1. Detweiler's system and its constraint amplification 



Detweiler's modification to ADM can be realized through one choice of the multipliers in (4.1) and (p~^). He 
found that with a particular combination the evolution of the energy norm of the constraints, + Ai^, can be 
negativ e definite when we apply the maximal slicing condi tion , K = 0. (We will comment more on his criteria in 



§IVC2.) His adjustment can be written in our notation in ( [4.l| ) and (4.2), as 

P„=-LaSj, (4.11) 

= La^iK,, - (l/3)if7,,), (4.12) 

= La^[3id(,a)S';) - {dia)j,jj''% (4.13) 

s^-',, =ic.3[5|=4-(l/3)7,,7'='], (4.14) 



everything else zero, where L is a constant. Detweiler's adjustment, ( |4.12| )-( 4.14 ), does not put constraint propagation 



equation to first order form, so we can not discuss hyperbolicity or the characteristic speed of the constraints. 

For the Minkowskii background spacetime, the adjusted constraint propagation equations with above choice of 
multiplier become 

'V'^MJ {-il/2)ih -iL/2)e5i - iL/6)hkJ y^^Mj ^ 
The eigenvalues of the Fourier transform are 



A' = A' = (-(i/2)fc2, -(L/2)fc^ -(4L/3)fc2 ± y'k^i-l + {4/9)L^k^ j). (4.16) 

This indicates negative real eigenvalues if we choose small positive L. 

We confirmed numerically, using perturbation on Minkowskii, that Detweiler's system presents better accuracy than 
the standard ADM, but only for small positive L. See the Appendix. 
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2. Differences with Detweiler's requirement 



We comment here on the differences between Detweiler's criteria for stable evolution and ours. 
Detweiler calculated the L2 norm of the constraints, Cp, over the 3-hypersurface and imposed the negative definite- 
ness of its evolution, 

Detweiler's criteria ^ '^t J CpC dV <0, V non zero Cp. (4-17) 

where CpC C^CpCa, and Gpa = d«ag[l,7y] for the pair of Cp = {TL,Mi). 

Assuming the constraint propagation to be dtCp = Ap'^Ca in the Fourier components, the time derivative of the 
L2 norm can be written as 

dtiCpCP) = {AP" + A'P + dtGP'')Cpd„. (4.18) 

Together with the fact that the L2 norm is preserved by Fourier transform, we can say, for the case of static background 
metric, 

Detweiler's criteria ^ eigenvalues of + are all negative V/c. (4.19) 

On the other hand. 

Our criteria <S4> eigenvalues of A are all negative V/c. (4.20) 
Therefore for the case of static background, Detweiler's criterion is stronger than ours. For example, the matrix 

A—(^^ where a is constant, (4-21) 

for the evolution system (Ci, C2) satisfies our criterion but not Detweiler's when \a\ > \/2- This matrix however gives 
asymptotical decay for (Ci, C2). Therefore we may say that Detweiler requires the monotonic decay of the constraints, 
while we assume only asymptotical decay. 

We remark that Detweiler's truncations on higher order terms in C-norm corresponds to our perturbational analysis; 
both are based on the idea that the deviations from constraint surface (the errors expressed non-zero constraint value) 
are initially small. 



D. Another possible adjustment 

1. Simplified Detweiler system 



Similar to Detweiler's (4.11), we next consider only the adjustment 

P,j =^ Koa-/,j, (4.22) 



all other multipliers zero in (4.1) and ( |4.2| ). 

On the Minkowskii background, the Fourier components of the constraint propagation equation can be written as 

[mj iv-(i/2)*fc, J [mj' ^^-^^^ 

and the eigenvalues of the coefficient matrix are 

A' = {0,0, Kok^ ± ,J k^{-l + Klk^)). (4.24) 

That is, the amplification factors become (0, 0, two negative reals) for the choice of relatively small negative kq. 
We also confirmed that this system works as desired. We give a numerical example in the Appendix. 
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2. Adjusting-Hamiltonian-constraints system 



Our final example is a combination of the one in § [VB| and that above, that is 

Pij = Koajij, (4.25) 
Rij = Kiajij, (4.26) 



all other multipliers zero in (4.1) and ( |4.2| ). Similar to the previous one, the Fourier transformed constraint propagation 
equation is 

which gives the eigenvalues 



^* J " U(l/2)»fc. - 2nA ' ) W)M, ' ^^'^^^ 



A' (0, 0, K^P ± y/k^{~l + Kok^ - 4ki)). (4.28) 

We can expect a similar asymptotical stable evolution by choosing kq and ki, so as to make the eigenvalues (0, 0, two 
negative reals). 



V. CONFORMAL-TRACELESS ADM SYSTEMS 



The so-called "conformally decoupled traceless ADM formulation" (CT-ADM) was first developed by the Kyoto 
group p5[. After the re-discovery that this formulation is more stable than the standard ADM by Baumgarte and 
Shapiro several groups began to use CT-ADM formulation for their numerical codes, and reported an advantage 
in stability . Along with this conformal decomposition, several hyperbolic formulations have also been proposed 

pOj-pH, but they have not yet been applied to numerical simulations. 

However, it is not yet clear why CT-ADM gives better stability than ADM. The Potsdam group |Q found that the 
eigenvalues of CT-ADM evolution equations has fewer "zero eigenvalues" than those of ADM, and they conjectured 
that the instability can be caused by "zero eigenvalues" that violate "gauge mode" . Miller applied von Neumann's 
stability analysis to the plane wave propagation, and reported that CT-ADM has a wider range of parameters that 
give us stable evolutions. These studies provide supports of CT-ADM in some sense, but on the other hand, it is also 
shown that an example of an ill-posed solution in CT-ADM (as well in ADM) 

Here, we apply our constraint propagation analysis to this CT-ADM system. 



A. CT-ADM equations 

Since one reported feature of CT-ADM is the use of the momentum constraint in RHS of the evolution equations 
[ p2t , we here present the set of CT-ADM equations carefully for such an replacement of the constraint terms. 

The widely used notation 16| is to use the variables {(f>,^ij ,K ,Aij ,r^) instead of the standard ADM {'^ij,Kij), 
where 

Itj ^er^'^ltj. (5.1) 
i,, = e-4^(i^,, - (l/3)7„if), (5.2) 
f - t)kl'\ (5.3) 
and we impose det7y = 1 during the evolutions. The set of evolution equations become 

(5.4) 
(5.5) 
(5.6) 







(-l/6)aif. 




{dt- 




-2aAy, 






- C0)K ^ 


a{l - Ki)R^^^ +a{l 




{dt- 




-e-4^(V,V,a)^^ + 
-fe-4^a(l/3)7,,K3[- 


e^^^aRf - e'^^a{l/3H,{l - K3)R'^^^ + a{KA,j - 



2Aii;A^ j 



(5.7) 



dtT' = -2{dja)A'^ - {4/3)K2a{djK)f^ + 12K2aA^'{dj(f>) - 2aAk^ {djf'') - 2K2ar''ijA^ kf' 
-2(1 - K2)a{d,Aki)f''f' + 2a(l - K2)A',P 

-dj {P'^dkf' - t'idkP') - i''\dkf3^) + {2/3)f^{dkl3'')) (5.8) 
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where Lp is the Lie derivative along the shift vector and i?^^^ is the 3-metric scalar curvature. Here we introduced 
parameters k which show where we replace the terms with constrai nts. F or e xample (aci,K2,K3) = (0, 0, 0) is the 
case of no replacement (the standard ADM equations expressed using (5.1)-(|5^), while Baumgarte-Shapiro |16| uses 
(ki,K2,K3) = (1, 1,0). 

The constraint equations in CT-ADM system can be expressed as 

H ^ e-^-^ijf^) _ 8e-'*'^f ^■(S.Sj-^) - Se-^-^f + 8e~*'^(a,0)r + (2/3)iC2 ^ A,^A'', (5.9) 

- (9,ifc.)7'-'' " (2/3)(a,A0 ~ + 6{d,^)Ah - r%Ah, (5.10) 

g'^r + dj^^\ (5.11) 

Here Ti.,M arc the Ha milt onian and momentum constraints and the third one, t/, is a consistency relation due to the 
algebraic definition of (^). 



B. Constraint propagation equations of CT-ADM 



Similar to th e ADM cases, we here show the propagation equations for (5.9)-( 5.11 ). The expressions are given using 



express R^j as 



( |3.5D and (|3.6[), but we have to be careful to keep using the new variable, F^, wherever it appears. Following ||16|, we 

R^^ = -(l/2)f "(az9„7„) + {l/2)%,d,r'' + (1/2)7^,9,^= + (l/2)f'=f(,,)fe 

+f "Tf,f,fe„ + 7'"Tf;.f,fc„ + 7'™f Lf fey . (5.12) 

The constraint propagation equations, then, are obtained by straightforward calculations as 

dtH = P^djH) - 2ae-'^'^f \d,Mj) + 2aKn - 2ae^'^'l'{d^f^)Mj - Aae-'^'>'{d,4))f^Mj - Ae-^'^f^{dja)M., 
+2K2e~^'''{d^a)f^Mj + 2n2e-^'^a{d^f^)Mj + 2n2e-^'^af^ {d.Mj) 

+16K2ae-'^'''{d^(t))f^Mj - {A/3)KiaKn, (5.13) 
^tM^ = -(l/2)a(ftH) + + aKM^ - {d^a)H - Ap'{^^4>)M, + p''j'\d^jik)Mj + id^pk)e-^^^''' Mj 

+ (1/3) (2ki + K3)(a,a)H + (1/3) (2ki + K3.)a{d,n) - 2K2aA^,Mj - il/3)K3ag^ ^jiH + 2K3a{d,(j))n, (5.14) 

dtG' = 2A'jg' + 2K2af''Mj. (5.15) 

These form a first order system, and the characteristic part can be extracted as 

/ n \ / /3' 2(-l + K2)aV^ o\ / n \ 

dtiMA -\ ((2/3)ki + (1/3)k3 - (l/2))a5^ /3'<5^ ] di [ Mj ] , (5.16) 

\g' J \ oj \g' J 

whose characteristic speeds are 

A' = (O, 0, 0, P\(3\P^± a^j'\l - K2)(l - (4/3)ki - (2/3)^3)) (no sum over /). (5.17) 

By analyzing the reality of the eigenvalues, the diagonalizability of the characteristi c ma trix, and the possibility of 
the symmetric characteristic matrix, we can classify the hyperbolic! ty of the system ( 5.16 ) as 

weakly hyperbolic <^ (1 - K2)(1 - (4/3)ki - (2/3)^3) > 0, (5.18) 
strongly hyperbolic <^ (1 - K2) = (1 - (4/3)ki - (2/3)^3) = 0, 

or (1 - K2)(l - (4/3)Ki - (2/3)k3) > 0, (5.19) 

symmetric hyperbolic ^ (-1 + K2) = (1 - (4/3)ki - (2/3)^3). (5.20) 

That is, for the non-adjusted system, (ki,K2,'^3) = (0,0,0), constraint propagation forms a strongly hyperbolic 
system, while the Baumgarte-Shapiro form gives only weakly hyperbolicity. (We note that the first-order version of 
CT-ADM by FrittcUi-Reula pO] has also well-posed constraint propagation equations. ) 
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C. Amplification factors on Minkowskii background 



For a Minkowskii background, the constraint propagation equations at the linear order become 



2(k2 - l)ikj 

((2/3)Ati + (1/3)^3 - (l/2))ifc, 

2k25'3 




The constraint ampHfication factor becomes 



A' = (0, 0, 0, 0, 0, ± ~ ^t2)(l - (4/3)ki - (2/3)k3)) 



(5.21) 



(5.22) 



That is, a' are either zero, pure imag inary or ± real numbers. For the non-adjusted system they are zero and pure 
imaginary (that is, the same as (3.11)), while the Baumgarte-Shapiro form gives us all zero eigenvalues. Therefore 
from our point of view, these two are not very different in their characterization of constraint propagation. 



VI. CONCLUDING REMARKS 



We have reviewed ADM systems from the point of view of adjustment of the dynamical equations by constraint 
terms. We have shown that characteristic speeds and amplification factors of the constraint propagation change due 
to their adjustments. We compared the equations for the ADM, adjusted ADM, conformal traceless ADM (CT- 
ADM) systems, and tried to find the system that is robust for violation of the constraints, which we can call an 
"asymptotically constrained" system. 

We conjectured that if the amplification factors (eigenvalues of the coefiicient matrix of the Fourier-transformed 
constraint propagation equations) are negative or pure-imaginary, then the system has better asymptotically con- 
strained features than a system they are not. According to our conjecture, the standard ADM system is expected 
to have better stability than the original ADM system (no growing mode in amplification factors). Detweiler's mod- 
ified ADM system, which is one particular choice of adjustment, definitely has good properties in that there are no 
growing modes in amplification factors. We also showed that this can be obtained by a simpler choice of adjustment 
multipliers. 

We also studied the CT-ADM system which is popular with numerical relativists nowadays. However, from our 
point of view, we do not see any particular advantages for CT-ADM system over the standard ADM system. 

The reader might ask why we can break the time-reversal invariant feature of the evolution equations by a particular 
choice of adjusting multipliers against the fact that the "Einstein equations" are time-reversal invariant. This question 
can be answered by the following. If we take a time-reversal transformation {dt —dt), the Hamiltonian constraint 
and the evolution equations of Kij keep their signatures, while the momentum constraints and the evolution equations 
of jij change their signatures. Therefore if we adjust 7ij-equations using Hamiltonian constraint and/or -equations 
using momentum constraints (supposing the multiplier has -I— parity), then we can break the time-reversal invariant 
feature of the "ADM equations". In fact, the examples we obtained all obey this rule. The CT-ADM formulation 
keeps its signature against the adjustments we made, so that we can not find any additional advantage from this 
analysis. 

Considering the constraint propagation equations is a kind of substitutional approach for numerical integrations of 
the dynamical equations. However, this might be one of the main directions for our future research, as Friedrich and 
Nagy pst impose the zero speed of the constraint propagation as the first principle when they considered the initial 
boundary value problem of the Einstein equations p^ . 

We are now applying our discussion to more general spacetimes, and trying to find guidelines for choosing ap- 
propriate gauge conditions from the analysis of the constraint propagation equations. These efforts will be reported 
elsewhere p7| . 
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APPENDIX A: NUMERICAL DEMONSTRATIONS OF ADJUSTED-ADM SYSTEMS 



We here show two numeric al dem onstrations of adjusted- ADM systems that were discussed in § IV C] (Detweiler's 
modified ADM system) and § [VD1 (s implified version). 

Detweiler's adjustment, (4.12)-(4.14), can be parametrized by a constant L, and our prediction from the amplifica- 
tion factor on Minkowskii background is that this system will be asymptotically constrained for small positive L. Fig.|l| 
is a demonstration of this system. We evolved Minkowskii spacetime numerically in a plane-symmetric spacetime, 
and added artificial error in the middle of the evolution. Our numerical integration uses the Brailovskaya scheme, 
which was described in detail in our previous paper The code passes convergence tests and the plots are for 401 
gridpoints in the range x — [0, 10], and we fix the time grid At ~ Q.2/S.X. The error was introduced as a pinpoint kick, 
in the form of A^j^j, = 10~^ at x = 5.0 and t = 0.25. We monitor how the L2 norm of the constraints [Ti.^ + Ai"^) 
behaves. From Fig^, we see that a small positive L reduces the L2 norm in time, which is the asymptotically con- 
strained feature we expected. The case of slightly larger L will make the system unstable. This is the same feature 
we have seen in the numerical demonstration of the A-system or adjusted-Maxwell/ Ashtekar systems p^ , for that 
case the upper bound of the multiplier can be explained by violation of the Courant-Friedrich-Lewy condition, while 
in this system we can not calculate the exact characteristics since the system is not first-order. 



Detweiler's adjustments 
on Minkowskii spacetime 




L= + 0.03 



1.0 
time 



FIG. 



IVC) 



1. Demonstration of the Detweiler's modified ADM system on Minkowskii background spacetime (the system of 
The L2 no rm o f the c onstraints is plotted in the function of time. Artificial error was added at t = 0.25. L is the 



parameter used in (4.12)-(4.14). We see the evolution is asymptotically constrained for small L > 0. 



Similarly, we plotted in Fig.gthe case of simplified version (the system of §[VD 1). We see the desired feature again 
by changing the parameter Kg that appear in (4.22). 
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Simplified Detweiler's adjustments 
on Minkowskii spacetime 




FIG. 2. Dem onstration of the simplified Detweiler's modified ADM system on Minkowskii back ground spacetime (the 
system of SIVD 1). For comparison with FigQ we set L — —kq, where kq is the parameter used in (4.22). We see the evolution 
is asymptotically constrained for small L > 0. 
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